ADA104293 


jJ 

.  J 


UNCLASSIFIED  1 

SECUH-TV  CLASSIFICATION  OF  THIS  PAGE  (Whon  Dolo  Enlormd)  \ 

u 

Cllpl 

X&l _ _ _ , 

|  )  .  REPORT  DOCUMENTATION  PAGE  1 

L 

Lit 

L  READ  INSTRUCTIONS 

BEKORE  COMPLETING  FORM 

1.  report  NUMBER - 

14092.20-M 

2.  GOVT  ACCESSION  NO 

ad-  A/o  f 

3  RECIPIENT'S  CATALOG  NUMBER 

293 

4.  TITLE  (md  SubtliU) 

\A  Boundary  Observability  and  Controllability 
for  Diffusion  Equations. 


7.  »UTHOH(.) 


'  Thomas  I.ISeidman 


/ 


».  PERFORMING  ORGANIZATION  NAME  AND  ADORESS 

University  of  Maryland 
Baltimore,  MD  21228 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 
Research  Triangle  Park,  NC  27709 


M.  MONITORING  AGENCY  NAME  »  ADOR£SS(lt  dllloront  Ina <  Controlling  Ollleo) 

/  - 


"S."  TYPE  OF  REPORT  a  PERIOD  COVERED 

/y  Final  Report » 

"  \tf  Dec  76  *  30  Jun  80  ,  ' 

6  pcrpowmimb  report  Mm«5eR 


B-  CONTRACT  OR  grant  NUMQERfftJ 

S' DAAG29-77-G  0061 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


12.  REPORT  0AT6 

Sp.p  fil 


U.  M OMBEft'OF  PAGES 

Ik _ 


IS.  SECURITY  CLASS,  (ol  Hit*  report; 

Unclassified 


ISo.  DECLASSIFICATION/' DOWNGRADING 
SCHEDULE 


la.  DISTRIBUTION  STATEMENT  (ol  Ihlm  Report; 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  tile  ebetr eel  entered  In  Slock  20,  II  dllloront  Iran  Report; 


NA 


§ 


DT 

ELEZ 

SEP  17  1981 


El  IT  *  '  T  7  r** 
w  t  ‘ 


18.  SUPPLEMENTARY  NOTES 

The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author (s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


IS.  KEY  WORDS  (Continue  on  rovormo  midm  It  nocooomry  mnd  identity  oy  oiock  numoor) 


|jfti  ABSTRACT  fCrttm  oxm—  Nfti  ft  nmcmemoey  mmd  Identity  ty  block  number)  — 

Investigation  has  focused  on  problems  of  observation  and  control  for  parabolic 
equations,  primarily  with  interaction  limited  to  the  boundary  of  the  region. 
Abstracts  and  bibliography  presented  herein  may  be  classified  under  six 
headings:  A.  Basic  theoretical  questions  on  control  of  parabolic  equations; 

U.  Problems  of  observation  and  identification  for  parabolic  equations;  C.  Othe 
system-theoretic  problems;  D.  Problems  and  considerations  (algorithms,  regularl 

?tc;)  related  to  approximation  and  numerical  computation-  E  111-r.^c  a  .  ,  I 
F.  froblems  involving  nnnlin.l.-  partial  differentia?  eau*fmncJ  P°sed  problemj 


DO 


FORM 
AM  71 


1473 


EDITION  OF  I  MOV  tt  IS  OBSOLETE 


UNCLASSIFIED 


T 


c 


ty, 

;  and 


<rrtimrv  r i  »wirir*Tini»  nr  tm»<  p»r.r 


Dafa  Fnfararfl 


«1  9  16  006 


// 


/V 


FINAL  REPORT 


ARO  grant  DAAG-29-77-G-006I 


Thomas  I.  Seidman  , 

Principal  Investigator 
University  of  Maryland  Baltimore  County 


"Boundary  observability  and  controllability  for  diffusion  equations" 


a.  Foreword 


b.  Problem  statements 


c.  Bibliographical  list 


d.  Summaries  of  results 


Accession  For 
trftS  GRAAI 
OTIC  TAB  “Cr 

Unannounced  □ 
Justification _ 

By - - 

Distribution/ 

Availability  Codes 
Avail  and/or 
1st  Special 


DTIC 

ELECTE 
SEP  17  1981 


DISTRIBUTION  J\i  Vi  W  •  .h" 

Approved  for  public  rdun-je; 
Distribution  Ui.'i>  i .t,.H 


a. 


Foreword  and  general  remarks 


This  grant,  as  extended,  covered  a  four  year  period  (essentially  1977- 1980). 
In  addition,  as  the  principal  investigator  was  out  of  the  country  (visiting  at 
the  Universite  de  Nice  for  the  academic  year  1980-  *81),  preparation  of  this 
final  report  has  been  somewhat  delayed;  since  much  of  the  work  of  this  additional 
half-year  is  based  on  and  is  a  continuation  of  work  done  during  the  grant  period, 
the  report  actually  includes  work  of  4%  years.  The  bibliographical  list  presents 
items  in  order  of  publication  or  expected  publication,  [l  ]  -  [34  ],  with  some 
brief  annotation.  The  summary  of  results  following  that  is  in  the  form  of 
abstracts,  prepared  for  this  report. 

During  the  grant  period  a  total  of  21  papers  and  reports  (items  [2  1  -  [23] 
noting  the  equivalence  of  r4  1 ,  [7])  were  prepared.  Of  these,  10  have  already 

appeared  or  have  been  accepted  for  publication  in  various  journals  while  8  others 
(The  3  appearing  in  the  annual  Seminaires  IRIA  have,  somewhat  arbitrarily,  been 
counted  with  these.)  are  in  conference  proceedings.  The  remaining  three  are  : 
a  note  [23 J  submitted  but  not  yet  accepted,  a  NASA  report  [8]  and  a  book 
chapter  [22]. 

Investigation  has  focused  on  problems  of  observation  and  control  for  para¬ 
bolic  equations,  primarily  with  interaction  limited  to  the  boundary  of  the 
region.  With  considerable  overlap,  the  items  [i]  -  [34]  may  be  classified  under 
six  headings: 

A.  Basic  theoretical  questions  on  control  of  parabolic  equations. 

[1 1,  K  ],  [5  J,  [13  ],  [14  J, [19  ],  [21  ],  [23  ],  [24  ],  [25  ],  [26  ],  [34  ]  . 

B.  Problems  of  observation  and  identification  for  parabolic  equations. 

[1],  [2  ],  [3  ]  [5  %  [6  ],  [11  ],  T18  ],  [22  ],  [27],  [33  ],  [34  ]. 

C.  Other  system-theoretic  problems. 

[2  ],[4],[20],  [32]. 

D.  Problems  and  considerations  (algorithms,  regularity,  etc.)  related  to  approx¬ 
imation  and  numerical  computation. 

[3], [5],  [6],  [11],  [12],  [13],  [14  ],  [16  ],  [18  ],  [19  ],  [22],  f26  ]. 

E.  Ill-posed  problems. 

[&  ]»  [8  1,  [9],  [11],  r  16],  [17],  [18  ],  [22  3,  [27  ]  . 

F.  Problems  involving  nonlinear  partial  differential  equations. 

[10  ],  [11  ],  [12  },  [14  ],  [15  ],  [24  y  [27  ],  [28  ],  [29  j,  [30  ],  [31  ],  f33  ],  p4  ]  . 


V 


As  general  observations,  note  the  following: 

(i)  The  original  proposal  focused  on  basic  theoretical  questions  of  observability 
and  controllability.  Much  of  the  published  work  is  relatively  abstract/ 
theoretical  in  style.  Nevertheless,  this  work  has  been  motivated  by  the  possi¬ 
bility  of  application  ( £ll 3  is  a  particular  case  in  point.)  and,  in  parti¬ 
cular,  the  heading  D  reflects  a  concern  with  considerations  relevant  to  success¬ 
ful  numerical  computation.  It  is  this,  of  course,  which. has  led  to  the  concerns 
with  regularity  of  optimal  controls  and  with  ill-posed  problems. 

(ii)  There  has  been  a  gradual  shift  toward  greater  concern  with  nonlinear 
problems.  In  part  this  reflects  an  attitude  that  the  'real'  problems  of 
application  are  nonlinear.  The  work  on  linear  problems,  of  course,  remains 
significant  -  not  only  as  a  model  but  because,  (e.g. ,  through  linearization) 
the  linear  results  apply  also  to  consideration  of  nonlinear  problems. 

(iii)  At  present  the  most  significant  outstanding  theoretical  problem  for  the 
linear  case  remains  patch  observability /controllability .  From  the  practical 
point  of  view,  the  outstanding  problems  are:  (1)  well-posedness  of  obser¬ 
vation  in  the  nonlinear  case  (for  which  this  is  no  longer  dual  to  exact 
null-control  ability),  (2)  improved  algorithms  for  computation  of  controls  and 
stabilizing  feedbacks  including  cases  with  constraints,  (3)  applications 
related  to  variational  inequalities  of  evolution  including  free  boundary 
problems,  (4)  realistic  modeling  of  stochastic  considerations. 


b.  Statement  of  some  problems. 

Consider  a  process  governed  by  the  parabolic  equation: 

(1)  u  +  Au  =  f 

where  A  is  a  second  order  elliptic  operator  —  say 

A  =  -  Au  or,  more  generally,  Au  =  -  7  •  a  Vu  +  qu 
with  specification  of  boundary  data  as 

(2)  uoru=n-avuorau+PuondQ. 

v  v 

If  (auv+  9°)  is  specified,  the  complementary  data  is 

(-  9u  +  au)/  (a2  +  B2  )  - 

v 

Equations  of  the  form  (1)  govern  a  wide  variety  of  phenomena  associated  with 
diffusion:  heat  conduction,  diffusion  of  a  solute,  etc. 

Distributed  control:  Select  f  so  as  to  minimize  a  '~ost'  of  the  form 

(3)  J  =  ||fl|2+  A  ||  u  -  u  H2+  ij.II  u  (T)  -  uj  /1 2 

(using,  say,  appropriate  l2  norms)  with  u  a  ’target  trajectory'  and  x 
a  'terminal  target'  given.  The  initial  and  boundary  data  are,  here,  also  given. 

Boundary  control:  Select  boundary  data  cp  so  as  to  minimize  (3)  with 

cp  for  f.  In  applications  this  is  a  typical  problem  as  one  often  has  no  possi¬ 
bility  of  interaction  except  at  the  boundary.  Indeed,  one's  interaction  mav  be 
limited  to  a  'patch'  on  the  boundary  (i.e.,  Ta  =  ’active'  boundary  c  dft  , 
with  the  constraint:  cp  =  0  on  d  Ta )  and  one  speaks  of  'patch  control'. 

Boundary  observation:  The  equation  and  boundary  conditions  are  assumed  known 
(e.g,,  homogeneous)  but  the  initial  state  unknown.  Given  observation  of 
the  complementary  boundary  data  as  well  (possibly  only  on  Ta  —  patch  observation 
one  wishes  to  use  this  to  determine  the  '•tate  u(T) . 

One  speaks  of  (boundary)  observability  if  this  last  is  possible.  It  is 
known  —  for  linear  problems  (the  corresponding  problem  for  nonlinear  equations 
is  an  interesting  target  for  future  investigation)  —  that  this  well-posed  if 
and  onl\  if  one  has  exact  nullcontrollability:  for  each  initial  state  there  is 
a  control  such  that  u(T)  =  0. 

In  these  contexts  one  seeks  information  on  the  existence  and  character¬ 
ization  of  optimal  controls,  continuity  of  dependence  on  various  aspects  of 
the  problem  and  methods  of  approximate  numerical  computation.  Major  results 
obtained  during  the  grant  period  (1)  extend  the  set  of  situations  for  which 
nullcontrollability  is  known  [1  ],  [  5  "]  , 

(2)  characterize  the  optimal  control  [5],  [19 ], [24  1  and  (3)  consider  aspects 
(regularity,  dependence,  etc.)  more-or-less  directly  related  to  computation 
[5 


Since  the  computational  aspects  lead  to  ill-posed  problems  —  e.g.,  the 
minimum  norm  boundary  control  from  u(0)  =  cu  to  u(T)  =»  u>  does  not  depend 
continuously  on  the  target  uS  —  it  has  been  relevant  to  consider  general 
computational  aspects  of  ill-posed  problems.  Such  problems  also  rise  in  connection 


with  system  identification,  e.g.,  determining  the  diffusion  coefficient  a 
in  specifying  A  above,  or  in  connection  with  problems  of  determining 
the  input  to  a  diffusion  system  from  observation  of  the  output  [  3 1, r  6], 
[  I6j,  [18]  .  The  Russian  school  (Tikhonov.,.)  has  been  particularly 
active  in  developing  computational  approaches  to  ill-posed  problems; 
these  have  been  of  great  importance  in  geophysical  contexts. 

Entirely  independent  of  the  above  is  a  nonlinear  diffusion-convec¬ 
tion-reaction  system  arising  in  semiconductor  theory  but,  with  modifi¬ 
cation,  applying  quite  gnerally  to  contexts  in  which  the  convection  is 
determined  by  the  electrostatic  field  which,  in  turn,  is  determined  by 
the  (unknown)  distributions  of  various  charged  species: 

uk  -  v  •  Jk  =  sk  (k  =  1.....K) 
where  the  'currents'  are  given  by 

Jk  =  \  7\  +  Uk\ 

with  'drift  velocities'  determined  by  the  field  E: 

Vk  =  3k  E  (qk  =  charge) 

and  E  =  -vt  for  a  potential  ^  given  by  the  usual  Poisson  equation: 

-A  ♦  =  N  +  (N=N(  •  )  given). 

Conservation  of  charge  gives  q.S  =  0. 

K  K  K 

A  discussion  of  the  interplay  between  theoretical  and  computational 
considerations  (for  the  steady  state  problem,  specifically  for  a  semi¬ 
conductor  device)  appears  in  fl2].  Major  results  [lO],  [30],  [31  ] 
are  on  existence  of  solutions  with  related  a  priori  estimates. 
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Summaries  of  principal  results  obtained; abstracts. 

The  principal  result  is  that  the  heat  equation  ut  = 

is  exactly  nullcontrollable  for  a  region  such  as 
indicated  with  boundary  control  restricted  to  the  patch 
1^  and  vanishing  on  SQ\ra  .  This  shows,  in  particular 
that  Russell's  'hyperbolic  =s  parabolic'  principle  is  far 
from  reversible  and  that  star-complementarity  is  far 
from  necessary  for  boundary  control  and  observation  of  the 
heat  equation.  For  a  general  patch  and  general  region  such  results 
unproved  but  state  prediction  from  such  observations  is  shown  to  be 
posed  given  an  a  priori  bound  on  initial  data. 


remain 
we  11- 


[2] 


Let  x(*)  satisfy  the  linear  equation 
observations  of  the  form  uu  =  2.^  x  ^ 


x  =  Ax  +  Bu 
(t^)  (finite 


and 

sum) 


consider  linear 
An  imprecise 


observation  vector  a;  is  given,  each  component  of  which  is  nominally  of 
this  form  (not  necessarily  linearly  independent),  and  one  seeks  an  opti¬ 
mum  (minimum  norm)  estimate  of  the  input  u( • )  .  The  equation  and  obser¬ 
vations  define  a  map:  (  F,u(*))  l — *  u>  :  (initial  data,  input)  I — *  (exact 
observation)  and  a  suitable  pseudo-inverse  is  to  be  constructed  relative 
to  an  indefinite  inner  product  which  ignores  §  . 


[3]  Let  u  be  known  to  satisfy  a  parabolic  equation  of  the  form: 
u  *  V  •  AVu  =  qu  in  (0,T)  x  W  with  au  +  £n’  Av u  =0  on 

(0,T)  X  and  initial  data  u  .  If  uo  and  the  coefficient  set 
o  =  (A,q,<j,  g)  are  (partially)  unknown,  one  asks  how  they  may  be 
determined  using  additional  boundary  observations.  (The  problems  are,  in 
general,  ill-posed  and  the  computational  considerations  are  further 
explored  in  [6  ],  [18  ',  [22  "]  .  ) 


[4  ^  For  general  linear  control  systems,  time-invariance  of  the  reachable 
set  is  a  purely  algebraic  corollary  of  exact  nullcontrollobility .  The 
result  is  formulated  for  time-dependent  problems:  Kg t  =  (states  at 
time  t  reachable  from  x  at  time  s)  is  independent  of  s,x  for 
s  s  r  <  t  if  every  state  at  r  is  controllable  to  0  at  t.  (This 
generalizes  a  result  of  Fattorini  and  was  used  by  Schmidt  in  proving  the 
bang-bang  principle  for  boundary  control  of  the  heat  equation.) 

[5]  Russell  had  shown  that  a  uniform  local  decay  rate  implies  boundary  con¬ 
trollability  for  the  wave  equation  u  =  A  u  and  that  this  implies 
corresponding  controllibility  for  u  =  Au  .  An  abstract  version  of  this 
is  proved  and  applied  to  give  controllability  for  a  class  of  variable 
coefficient  hyperbolic  and  parabolic  problems.  The  principal  new  results 
for  the  parabolic  case  are:  (1)  a  characterization  of  minimum  norm  null- 
controls  as  complementary  boundary  data  of  the  adjoint  equation  and  (2) 
continuous  dependence  of  the  optimal  control  (in  Hs((0,T')  X  dft  for 
0  <T/  <T  and  weak  in  L2((0,T)ydQ))  on  the  coefficients  of  the  equa¬ 
tion.  (The  first  of  these  gives  analyticity  in  t  of  optimal  nullcon- 
trols;  compare  [19],  ) 


T’.fc*.. 


A  number  of  problems  are  discussed  briefly.  Principal  new  results  are: 

(1)  an  abstract  result  on  computing  solutions  of  linear  problems  with 
convex  constraints  applied  to  optimal  boundary  control  and  (2)  uniqueness 
for  some  problems  of  identifying  q  in  u  =  u  -  qu  . 

See  [4  ]  . 

There  is  a  considerable  variety  of  approaches  to  ill-posed  problems 
(as  geophysical  inverse  problems)  which  have  significant  differences  in 
method  and  interpretation  based  on  varying  assumptions  regarding,  e.g., 
the  nature  of  the  measurement  uncertainties.  This  primarily  expository 
report  compares  various  approaches  to  a  posteriori  interpretation  of 
imprecise,  inadequate  data  and  gives  uncertainty  estimates  where  possible. 

Let  Ax  =  b  be  a  linear  problem  (here  presented  abstractly).  The 
method  of  'least  squares'  is  often  applied  by  writing  x  =  72  CL^  e, 

(|e  [  an  orthonormal  basis),  taking  enough  terms  k  =  1,  . . . ,K 

to  permit  good  approximation  and  then  determining  x^=  J ^  to  mini¬ 

mize  ||Ax^-bjP.  For  an  ill-posed  problem  (A  ^  unbounded)  examples  are 

constructed  for  which  an  exact  solution  exists,  exact  computation  is  assumed 
but  jx^  |  is  unbounded  as  K  is  increased.  Contrary  to  the  situation 

for  well-posed  problems,  the  approximation  can  get  worse  as  more  terms 
are  taken. 

A  general  reaction-diffusion-convection  system  for  charged  species  is: 

-  V -DkCvuk+ qk  ufc  Vtlr]  =  Sk  (k=l,...,K 

-  A  *  -  N  +  ^  q  uk 

where  q  is  the  charge  and  u,  tne  concentration  of  the  k-th  species  and 

K  K 

t|r  is  the  electrostatic  potential.  Under  suitable  hypotheses,  existence  of 
solutions  (with  each  u^  positive)  can  be  proved. 

In  particular,  this  applies  to  a  specific  system  of  semiconductor  theory  with 
D  dependent  on  |  |  and  S  =  -R  given  by  the  Shockley-Reed-Hall 

K  K 

model  of  recombination. 

Let  <p=  cp  (t)  be  a  thermal  transient  of  interest  and  suppose  one  can 
observe  the  output  ijr(t)  =  u(0,t)  determined  by  the  diffusion  equation 
u  =  [Du  ]  with  u(l,t)  =  cp(t),  u  (0,t)  =  0  and  some  initial  con- 

t  XX  X 

dition  u  =«  u^.  For  D  =  D(x)  or  for  D  =  D(u)  analytic  in  u  it  is 

shown  that  cp  is  uniquely  determined  by  exact  knowledge  of  $  .  Some 
computational  algorithms  for  determining  <p  from  \|i  are  proposed 
and  convergence  results  stated.  (This  is  related  to  the  nroblem  of  deter¬ 
mining  temperature  transients  in  a  gun  at  firing;  there  D  =  D(x,u), 
after  taking  polar  coordinates.) 

The  relation  is  given  between  the  existence  result  of  TO  1  and  an 
earlier  paper  (1972)  on  numerical  computation.  Since  the  result  uses  a 
substitution  vk  =  uk  exp[qk^],  an  explicit  bound  on  uk  in  terms  of 


the  voltage  across  the  device  requires  estimating  the  solution  of 
-  A  <P  =  Xe  as  X  increases  (compare  [l5[).  A  calculation  is  also 
presented  relevant  to  the  question  of  bifurcation  as  voltage  increases. 

[13  2  It  is  shown  that  an  infinite  horizon  optimal  control  problem  has  a  re¬ 
formulation  as  construction  of  a  stabilizing  feedback.  In  this  form  one 
has  a  pure  initial  value  problem  giving  a  semigroup  which  is  shown  to  be 
holomorphic. 

[14 2  After  a  summary  of  the  argument  of  [19  [,  q.v.,  asymptotic  estimates 
are  obtained  as  t  -»  0  +,T  -  for  the  H  (  gQ)-  norm  of  cp(t)  where  ^ 
is  an  optimal  Neumann  boundary  control  for  the  heat  equation  on 
(0,T)xf2.  Similar  results  are  obtained  for  a  semilinear  equation.  (These 
estimates  are  related  to  the  convergence  rates  obtainable  in  various 
methods  of  numerical  computation.) 


[15  2  Generalizing  an  argument  of  [12  ]  for  f(u)  =  e 
f=e"h  with  h  and  convex,  one  has  max  |u|  = 

as  x  00  •  (For  e'u  this  gives  log  X  growth, 
discovered  subsequently,  shows  that  is  sharp.) 


it  is  shown  that,  for 
(T(f'1(l/  X)) 

An  explicit  solution. 


[16  2  A  comparison  is  given  of  numerical  methods  for  computing  optimal  bound¬ 
ary  controls  for  the  heat  equation  in  three  contexts:  minimum  norm  controls 
with  no  constraint  or  nonnegativity  constraint  (semi -const rained:  (V.ad 
unbounded)  and  time-optimal  controls  with  a  given  bound.  The  method 
recommended  in  general  was  a  projection  method  (FEM)  with  regularization  . 

An  argument  was  given  to  show  this  projection  method  converged  even  with¬ 
out  regularization  when  applied  to  unconstrained  nullcontrol.  (Note: 

The  comparison  did  not  include  switching-point  methods  applicable,  noting 
the  bang-bang  principle,  to  certain  time-optimality  problems  in  one  space 
dimension. ) 


[l7},[18]  These  were  prepared  originally  in  1978,  based  on  still  earlier  work, 
and  have  been  extensively  distributed  in  report  form.  The  principal 
theoretical  results  [17 3  are:  (1)  a  general  convergence  theorem  for 
approximations  F^(x)«  b^  to  F(x)  =  b  .  (The  correct  notion  of  conver¬ 
gence  F^-»  F  is  'graph  subconvergence'  in  x  Y  :  if  In 

graph  (F-k)  and  x^-^-x,  y^_>.y,  then  (x,y)  in  graph  (F).)  and  (2) 

discussion,  with  error  estimates,  of  the  method  of  'generalized  inter¬ 
polation',  particularly  in  the  linear  case  in  Hilbert  space  (approxima¬ 
tion  x^  to  Ax  =  b  given  by  <  x^,  A  j) j  >  =  <b,T]j>  for 

j  =  1, . . .  ,K  with  Xj^  in  sp  JA  *7]  ^  )  . 


The  applications  considered  [18  ]  are  (1)  the  inversion  of  Nemytsky 
operators (nonlinear  back-substitution),  (2)  recovery  of  a  diffused  signal 
(compare  written  later)  and  (3)  system  identification  for  systems 

governed  by  parabolic  and  functional  differential  equations  (compare  [3], 
written  earlier). 


or  u 


[19] 


Let 


ut  +  Au  = 


'an 


data  and  suppose 

+  'j.  ||ucr)-to|^ 


0  ( e . g .  ,  A= -  £ )  , 

Cp  is  optimal,  minimizing  J 

with  the  target  trajectory  u 
.  s. 


=  <P 
=  I! 


with  given  initial 


W  E 

analytic  in 


+  A||U-  u 


t. 


H2 

Q 

Then 


is  analytic  in  t  as  an  H  (  gQ)  -  valued  function  (s  as  large  as  is 
consistent  with  the  regularity  of  0  ft  ,  d  and  the  coefficients  of  A) 
on  a  complex  domain  including  the  open  real  interval  (0,T).  This  is  also 


true  for 


(exact  control  to  cu  assuming  this  possible)  and  for 


infinite  horizon  (stabilization)  problems;  compare  [5],  [13  ],  [14]. 


[20  J  A  variety  of  situations  can  be  modeled  by  this  problem  or  more  compli¬ 
cated  variants:  one  seeks  to  maximize  the  expected  discounted  return 
associated  with  a  choice  between  a  'gamble*  with  known  success  probability  p 
and  another  with  (unknown)  probability  p.  Beginning  with  an  assumed  0 
distribution  for  p,  one  plays  this  so  long  as  it  provides  the  greater 
expectation  (with  Bayesian  update  of  the  distribution).  This  is  thus  an 
'optimal  stopping  time'  problem  with  optimal  strategy  characterized  by 
the  'free  boundary'  separating  the  choice  regions  in  parameter  space. 

An  approximation  is  obtained  for  the  free  boundary. 


[211  Consider,  e.g.,  the  trace  rju  =  u  for  eigenfunctions:  -  4  u  *  W  in  ft , 
uv=  0  on  30  with  |lu||^=  1.  It  is  shown  that  ||fi>|)p  is  bounded  below 

if  F  is  a  large  enough  part  of  30.  For  <u  =  uv|^  and  -  ^u  =  Au,  u  |  •  0 

one  has  |)u>  ]|  £  )  •  The  argument  is  an  application  of  known  results 

in  boundary  control  theory  for  utt  =  Au .  Other  related  results  are  also 
obtained. 


[22]  This  is  primarily  expository,  addressed  to  an  audience  of  'engineers'. 
The  discussion  of  computational  approaches  to  problems  of  observation, 
control  and  identification  for  parabolic  equations  focuses  on  the  ill-posed 
problems  which  may  arise  and  the  implications  of  this  for  computation. 


[23] 


One  seeks  to  minimize  J  = 


IcpUj,  +  Alfu  (T)-£iln 


where 


ut=  Au  +  f(u)on 


Q  =  (0  ,T)x  Q ,  u  =  cp  on  £=  (0,T)  y  3Q  and  u(0)  =*  <u  .  Note  that,  even 
for  f  ;0,  J(  q,)  need  not  be  finite  for  cp6L2(v)  and  J  is  not  lower 
semi-continuous.  Under  suitable  hypotheses  on  f,  a  technical  trick  — 
apparently  applicable  also  to  a  variety  of  other  problems  —  shows  that 
if  j(  (p^  ,  u^) }  is  a  minimizing  sequence  with  (  .u^)  then 


(d^.u*)  minimizes  J  admissibly  (including  the  consideration  of  constraints, 
not  necessarily  convex  with  respect  to  u  )  . 


[24]  Let  u^  +  Au  +  f(u)  =  cp  with  u(0)  =  0,  u|  =  0  and  ©  to  minimize 

J  =  [f©J^  +  \||u- u||2  +  p, |f u  (T)-<i||2  .  Assume  f’=g  with  Osgi  K.  Then  (1) 

optimal  controls  exist,  (2)  cp  I — *  u  is  Frechet  differentiable  (This  is 
shown  here  only  under  dimension  limitations:  OclRm  with  m  <  10  for 

A  =  -A  *  but  has  since  been  shown  to  hold  in  general.)  and  (3)  if  the 


2 

data  u,u>  is  'small  enough'  and  ^  <  X  then  the  optimal  (cp*,u*)  is 
unique  and  depends  continuously  on  u,v:  and,  in  a  certain  sense,  on  f. 
Boundary  control  is  also  considered. 

[ 25  i  The  first  result  is  exact  nullcontrollability  by  boundary  control  for 

certain  diffusion  equations  in  one  space  dimension  with  time-dependent 
coefficients;  e.g.,  for  u  =  u  -  q(t,x)u  on  (0,T)x  (0,1)  provided 

q  is  analytic  in  t  (at  t  =  0)  uniformly  in  x£(0,l).  The  second 
result  is  that  when  a  parabolic  equation  is  shown  to  be  exactly  boundary 
nullcontrollable  (so  there  exists  Initial  data  optimal  control  to 
give  u(T)  =  0),  then  log  )|<Xp |j  =  ^(l/f)  as  T  — - ♦  0  +  ;  This  shows  how 

much  harder  one  must  work  to  control  quickly.  Analogous  results  hold 
for  the  dual  observation  problems. 

£26  The  principal  result  addresses  the  following  difficulty  which  arises, 
e.g.,  in  justification  of  optimality  conditions  for  control  problems. 

Given  a  minimum  x  of  J :  X  -♦  R  one  may  be  able  to  justify  differen¬ 
tiability  so  j‘  (x)  =  0  if  one  knows  xg  YCX  (regularity)  while 
the  condition  J'(x)  =  0,  obtained  formally,  may  imply  x£Y  —  e.g., 
if  j‘  has  the  form  y  -  K(y)  with  K:X->  Y.  It  is  shown  that  this 
circular  argument  may  be  rigorously  justifiable,  using  Ekeland's 
approximate  variational  principle. 

[27  1  Various  nonlinear  diffusion  processes  are  governed  by  equations  of  the 
general  form  ufc  =  v*k(u)  (Examples  of  such  problems  are  seepage 

of  water  or  oil  underground  and  heat  conduction  with  a  temperature-depend¬ 
ent  coefficient.)  To  identify  the  particular  nonlinearity  k(>)  from 
boundary  observations  is,  in  general,  an  ill-posed  problem.  For  0  <  k 
analytic  and  Dirichlet  data  <p=  cp(x  )  (smooth  but  not  constant), 
observe  the  resulting  flux  %  =  k(u)  yu  -  n  on  (0,T)x  dQ.  It  is  shown 
that  the  pair  (  cp,  )  uniquely  determines  k(’).  (It  is  anticipated  that 
a  similar  result  holds  if  the  flux  is  specified  and  the  trace  observed. 
Some  specific  numerical  procedures  may  be  included  as  well.) 

[28  T  If  f:  1R  -*  IR  is  bounded  above  as  r  -*  <*■  and  below  as  r  _*  -  ®  and 
if  A  provides  a  suitable  maximum  principle  (e.g.,  A  =  -  a  with 
bounded  Dirichlet  data),  then  Au  =  f(u)  has  a  solution. 

[29  J  Consider  A:u  i — v  .  a(-,|vu|)yu  with  ra(-,r)  having  suitable 

growth  properties.  Coercivity  estimates  are  obtained  giving  existence, 
uniqueness  and  continuous  dependence  results  for  solutions  of  related 
elliptic  and  parabolic  problems.  An  example  would  be 

a(-,r)  =  r  ®/[<u(  *  )  +  ra  ]  with  a  >  0,  Ojiu  6L^.  An  application  given 
is  to  generation  of  periodic  eddy  currents  in  an  inhomogeneous  nonlinearlv 
magnetic  material. 


[30  ^  pi  "*  The  system  governing  diffusion  and  convection  of  hole  and  electron 
concentrations ,  with  Shock  ley -Reed -ha  11  recombination,  in  a  semiconductor 
device  (compare  [10]  for  the  steady  state  problem)  is: 

ut-  v  -aC  v 11  +  «  v  i1  =  R  , 

v  -  V  ■  b  [  Vv  -  v  7  if  1=  -K, 

&l|r=N  +  U'V, 

2 

R  =  [  u  v  -  n  1/  t  ,  i  =  t  (  ■  ,u ,  v)  . 

Existence  results  are  obtained  for  positive  solutions  of  this  system.  The 
argument  of  [30'’  is  for  a,b,T  constant  ,  n^  =  1.  The  argument  of  [31] 
considers  a,b  functions  of  x,  j  and  r=  t  (  •  ,  ts ,  v,  |7  ^1),  although 
under  somewhat  different  hypotheses  than  in  [10  It  is  anticipated  that 
this  can  be  coupled  with  heat  conduction  and  permits  n  =  n(x),  not  constant. 

[32  '  For  C  closed  and  convex  in  a  Hilbert  space  H  and  S  a  subspace 

one  considers  when  the  nearest  point  to  0  in  CHS  is  also  the  nearest 
point  in  C  to  some  x  £  S -i-.  This  is  related  to  a  representation  theorem 
for  certain  optimal  control  problems  with  constraints  as,  e.g.,  nonnegative 
control. 

p3  ]  The  determination  of  the  state  for  a  Stefan  problem  from  observation 
of  the  motion  of  the  free  boundary  is  a  well-pooed  problem. 

p4  ]  A  book  in  three  main  parts:  (l)  Elliptic  Equations,  (2)  Parabolic 

Equations,  (3)  Observation  and  Control  of  Parabolic  Equations.  The  empha- 
sic  in  (i)  is  on  variational  methods  with  applicability  to  nonlinear 
problems.  The  emphasis  in  (2)  is  on  semigroup  methods  and  semi  linear  equa¬ 
tions  with  inhomogeneous  boundary  conditions.  Pari  (3)  treats  various 
topics,  concentrating  primarily  on  boundary  observation  and  control  along 
the  lines  of  [5  J,  [13  ],  f  19  } ,  [23],  [24],  etc. 


